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Abstract 

An  elastic-viscoplastic  model  for  describing  the  anisotropic  high-strain  rate  behavior  of  both 
low-symmetry  and  high-symmetry  textured  materials  is  proposed.  Yielding  is  described  using  a 
recently  developed  criterion  which  can  capture  simultaneously  anisotropy  and  compression- 
tension  asymmetry  associated  with  deformation  twinning.  The  anisotropy  coefficients  as  well  as 
the  size  of  the  elastic  domain  are  considered  to  be  functions  of  the  accumulated  plastic  strain. 
The  specific  expressions  for  the  evolution  laws  are  determined  using  a  multi-scale  methodology, 
i.e.  experimental  measurements  of  crystallographic  texture  and  uniaxial  stress-strain  curves,  poly¬ 
crystalline  calculations,  and  macroscopic  scale  interpolation  techniques.  An  overstress  approach  is 
used  to  incorporate  rate  effects  in  the  formulation.  Applications  of  the  model  to  the  description  of 
the  high  strain-rate  response  of  low-symmetry  (hexagonal-close-packed  zirconium)  and  high-sym- 
metry  (body-centered-cubic  tantalum)  pre-textured  metals  are  presented.  The  very  good  agreement 
between  the  simulated  and  experimental  post-test  geometries  of  the  Taylor  impact  specimens  in 
terms  of  major  and  minor  side  profiles  and  impact-interface  footprints  shows  the  ability  of  the 
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model  to  describe  with  fidelity  the  differences  in  the  evolution  of  anisotropy  between  zirconium 
and  tantalum. 

©  2006  Elsevier  Ltd.  All  rights  reserved. 
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1.  Introduction 

The  importance  of  an  adequate  description  of  plastic  anisotropy  has  been  demon¬ 
strated  in  many  low  strain-rate  forming  applications  (e.g.  Cazacu  and  Barlat,  2003;  Yoon 
et  al.,  2006;  Soare  et  al.,  in  press,  etc.).  Recent  efforts  in  the  development  of  computational 
models  to  describe  temperature  and  strain-rate  effects  on  the  inelastic  response  have  given 
rise  to  robust  predictive  methods  for  simulating  the  anisotropic  high  strain-rate  behavior 
of  body-centered-cubic  (bcc)  and  face-centered-cubic  (fee)  poly  crystals.  Maudlin  et  al. 
(1999a,b)  employed  an  anisotropic  representation  of  the  yield  surface  obtained  from  poly¬ 
crystal  calculations  to  model  the  mechanical  response  of  a  pre-textured  bcc  tantalum 
sheet.  Anand  et  al.  (1997)  developed  finite-element  (FE)  poly  crystal  models  for  describing 
the  high  strain-rate  and  large  deformation  behavior  of  both  bcc  and  fee  polycrystals. 
These  models  describe  qualitatively  well  the  important  mushrooming  of  the  Taylor  impact 
specimens.  In  contrast,  the  prediction  of  the  geometry  changes  observed  in  hexagonal- 
close-packed  (hep)  specimens  remains  a  challenge  (Maudlin  et  al.,  1999b).  Metallographic 
investigations  have  shown  that  while  deformation  in  hep  materials  under  quasi-static  and 
dynamic  loading  includes  both  slip  and  deformation  twinning,  the  propensity  of  twinning 
is  higher  for  high-rate  loading  and/or  low-temperature  deformation.  Therefore,  in  order 
to  describe  with  accuracy  the  dynamic  response  of  hep  materials  it  is  imperative  to 
account  for  the  most  important  sources  of  anisotropy  in  the  given  material:  slip  and/or 
twinning  activity,  substructure  evolution  at  grain  level,  and  texture  development  during 
deformation. 

Since  in  crystal  plasticity  models  the  distribution  of  crystal  orientations,  the  available 
slip  and/or  twinning  deformation  systems,  and  the  stress  levels  necessary  to  activate  them 
are  taken  into  account  explicitly,  the  evolution  of  anisotropy  due  to  texture  development 
can  be  characterized  by  measuring  the  initial  texture  and  calculating  grain  reorientation 
(i.e.  updating  the  texture)  using  a  suitable  homogenization  scheme.  Recently,  the  develop¬ 
ment  of  crystal  plasticity  models  for  hep  metals  and  the  incorporation  of  crystal  plasticity 
calculations  directly  into  finite  element  (FE)  codes  have  received  much  attention.  Models 
that  account  for  both  slip  and  twinning  activity  and  employ  Taylor  (Bishop  and  Hill, 
1951)  or  self-consistent  averaging  schemes  (Lebensohn  and  Tome,  1993)  to  predict  the 
aggregate  behavior  have  been  proposed  (e.g.  Staroselsky  and  Anand,  2003;  Salem  et  al., 
2005;  Tome  et  al.,  2001;  Tome  and  Lebensohn,  2004  etc.).  So  far,  these  models  have  been 
validated  only  for  quasi-static  deformation  processes. 

This  paper  presents  a  macroscopic  elastic- viscoplastic  model  for  describing  the  dynamic 
anisotropic  plastic  response  of  textured  metals.  The  aim  is  to  develop  a  model  that  can  be 
applicable  to  both  low-symmetry  and  high-symmetry  pre-textured  materials.  Key  in  this 
development  is  the  use  of  a  recently  proposed  anisotropic  yield  function  (Cazacu  et  al., 
2006).  This  yield  function  (denoted  in  what  follows  as  CPB06)  is  capable  of  describing 
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simultaneously  anisotropy  and  tension/compression  asymmetry.  Moreover,  if  the  param¬ 
eter  associated  with  strength- differential  effects  is  set  to  zero,  the  CPB06  criterion  also 
describes  accurately  the  yielding  of  metals  with  no  tension/compression  asymmetry. 

The  plan  of  the  paper  is  as  follows.  Section  2  presents  an  overview  of  the  proposed  elas¬ 
tic- viscoplastic  model.  In  Section  3,  numerical  aspects  related  to  the  FE  implementation  of 
the  model  are  discussed.  In  Section  4  we  present  an  application  of  the  proposed  formula¬ 
tion  to  hep  materials.  In  particular,  we  describe  how  the  parameters  involved  in  the 
evolution  laws  for  a  high-purity  zirconium  clock-rolled  plate  can  be  determined,  based 
on  experimental  information  (texture,  quasi-static  uniaxial  tension  and  compression  load¬ 
ing  curves  and  metallographic  evidence)  and  simulation  results  obtained  with  the  visco¬ 
plastic  self-consistent  (VPSC)  poly  crystalline  model  of  Lebensohn  and  Tome  (1993). 
Next,  a  validation  of  the  anisotropic  low  strain-rate  formulation  is  provided  by  comparing 
FE  simulations  with  experimental  results  on  the  zirconium  plate.  The  elastic- viscoplastic 
extension  of  the  model  is  then  used  to  simulate  Taylor  cylinder  impact  test  results  for 
the  same  material.  We  conclude  with  an  application  of  the  model  to  bee  materials  perform¬ 
ing  the  analysis  of  the  Taylor  cylinder  impact  test  on  pre-textured  tantalum.  It  is  shown 
that  the  model  describes  with  accuracy  the  striking  difference  between  the  geometries  of 
the  post-test  specimens  of  zirconium  and  tantalum. 

2.  Anisotropic  elastic-viscoplastic  model 

Our  objective  is  to  develop  a  macroscopic  anisotropic  model  that  describes  simulta¬ 
neously  the  influence  of  strain-rate,  temperature,  and  evolving  texture  on  the  inelastic 
response  of  textured  metals  for  monotonic  loading  paths.  Thus,  an  elastic/viscoplastic 
approach  (see  Perzyna,  1966)  is  adopted.  The  basic  assumption  is  that  the  viscous  proper¬ 
ties  of  materials  become  manifest  only  after  the  passage  to  the  plastic  state.  Thus,  the 
strain-rate  g  can  be  decomposed  additively  into  an  elastic  part  (gE)  and  a  viscoplastic  part 
(gvp),  the  latter  representing  the  combined  viscous  and  plastic  effects,  i.e. 

g  —  T  ^vp  ( 1 ) 

The  evolution  of  the  viscoplastic  strain-rate  is  considered  to  be  given  by  an  overstress 
type  law  of  the  form: 


evp  =  y 


(2) 


where 


for  /  ^  0, 
for  /  >  0. 


(3) 


In  Eq.  (2), /is  the  quasi-static  yield  function,  g  is  the  quasi-static  plastic  potential,  y  a 
viscosity  parameter,  and  m  is  a  strain-rate  sensitivity  constant. 

The  overstress  law  (2)  reduces  to  the  classic  plastic  flow  rule  if  the  strain  rate  goes  to 
zero  (Perzyna,  1966).  The  quasi-static  yield  function  is  the  CPB05  yield  function  (Cazacu 
et  al.,  2006): 


(4) 


1004 


B.  Plunkett  et  al.  /  International  Journal  of  Plasticity  23  (2007)  1001-1021 


where  Y  is  a  hardening  parameter,  which  is  a  function  of  the  accumulated  equivalent 
viscoplastic  strain  svp,  and,  possibly,  of  the  temperature  T,  while  o  is  an  equivalent  stress 
defined  in  the  CPB05  framework  as  (Cazacu  et  al,  2006;  Plunkett  et  al.,  2006): 

o  =  B[(|^  |  -  kZx)a  +  (|272|  -  kZ2)a  +  (l^l  -  kZ3)a]'  (5) 

with  B  being  a  constant  defined  such  that  o  reduces  to  the  tensile  yield  stress  in  the  rolling 
direction,  i.e. 


L(I*1 1  -  W  +  (\02\  -  k$2)a  +  (|<p3|  -  k03 )fl 

and 


In  Eq.  (5),  Zx,  Z2,  are  the  principal  values  of  the  tensor  Z  =  L [S],  where  L  is  a  fourth- 
order  orthotropic  tensor  that  accounts  for  the  plastic  anisotropy  of  the  material  and  S  is 
the  deviator  of  the  Cauchy  stress  tensor,  A;  is  a  strength-differential  parameter,  and  a  is  a 
homogeneity  constant.  Note  that  for  the  case  of  a  =  2,  k  =  0,  and  L  equal  to  the  fourth- 
order  identity  tensor,  CPB05  reduces  to  the  von  Mises  criterion  and  B  reduces  to  ^/3/2, 
i.e.  the  constant  associated  with  the  von  Mises  effective  stress.  The  effective  plastic  strain 
associated  with  this  anisotropic  yield  function  is  calculated  using  the  principle  of  equiva¬ 
lence  of  plastic  work  (Hill,  1987). 

In  the  above  equations  all  stresses  and  strains  are  expressed  in  a  materially  embedded 
coordinate  system  such  as  to  ensure  that  the  resulting  incremental  constitutive  law  is 
objective  (Yoon  et  al.,  1999).  The  evolution  of  anisotropy  due  to  evolving  texture  is 
accounted  for  by  means  of  a  multi-scale  procedure,  which  essentially  consists  in  using 
experimental  textures  and  stress-strain  curves,  polycrystalline  calculations,  and  macro¬ 
scopic  scale  interpolation  techniques.  The  procedure  starts  by  fitting  the  parameters  of  a 
polycrystal  model  (we  use  the  VPSC  formulation)  to  reproduce  the  mechanical  response 
of  the  pre-textured  material  along  a  given  deformation  path  (note  that  the  measured  initial 
texture  is  needed  as  input  of  the  polycrystal  model).  Specifically,  using  VPSC,  the  poly¬ 
crystal  is  pre-strained  up  to  a  given  deformation  level  W  along  a  given  deformation  path. 
Then,  in  order  to  quantify  the  induced  anisotropy,  the  pre-strained  polycrystal  is  further 
“numerically”  probed  along  various  directions.  The  calculated  yield  stresses  along  these 
various  test  directions  are  in  turn  used  to  determine  the  coefficients  L (e(p), A(e(p), a(g(p) 
involved  in  the  macroscopic  CPB05  yield  criterion  and  thus  to  obtain  the  macroscopic 
yield  surface  corresponding  to  the  texture  for  the  given  level  of  pre-strain,  by  calculating 
aj  according  to  Eq.  (5),  as  well  as  Yj  =  7(e(p).  This  procedure  is  repeated  for  a  finite  set  of 
pre-strain  levels,  say  slyp  <  <. . .  <  syp.  Further,  an  interpolation  procedure  is  used  to 

obtain  the  macroscopic  yield  surfaces  corresponding  to  any  pre-strain  level  (Plunkett 
et  al.,  2006),  i.e. 

ff(<T,eVp)  =  ^(evpK  +  [1  -  £(eVp)]ff;'+1 

F(£VP,  T)  =  £(£vp)  •  YJ  +  (1  -  £(evp))  •  F+1 


(8a) 

(8b) 
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for  any  e(p  evp  ^  e£‘,  j  =  1 . . .  n 
Eq.  (8)  is  defined  as 


£(fivP) 


pj+ 1 


Gy 


i+ 1 


£yp 


1.  The  interpolation  parameter  £(gvp)  appearing  in 


(9) 


It  is  worth  noting  that,  although  the  above  procedure  delivers  a  viscoplastic  response 
for  any  applied  state  at  any  accumulated  deformation,  such  mechanical  response  repre¬ 
sents  a  good  approximation  of  the  actual  behavior  of  the  material  only  for  applied  load¬ 
ings  with  moderate  deviations  with  respect  to  the  adopted  pre-strain  deformation  path. 
For  quasi-static  strain  rates,  the  proposed  elastic/viscoplastic  model  reduces  to  an  elas¬ 
tic/plastic  model  (see  Perzyna,  1966). 


3.  Integration  algorithm 

In  order  to  simulate  high-rate  deformation  processes,  the  elastic-viscoplastic  model 
described  by  Eqs.  (1)— (7),  (8a),  (8b)  was  implemented  in  the  commercial  FE  code  ABA- 
QUS  (2003).  The  algorithmic  aspects  related  to  the  FE  implementation  are  presented  in 
what  follows. 

In  displacement-based  FE  formulations,  stress  updates  take  place  at  the  Gauss  points 
for  a  prescribed  nodal  displacement.  We  start  from  the  known  converged  state 
[ff ,  £yP5  •>  ^yp]  and  calculate  the  corresponding  values  at  t  +  At:  [^+A?,  £^A',  ot+At,  e^].  In 
this  incremental  process,  the  total  strain  increment  As  is  decomposed  into  an  elastic  part 
A eE  part  and  a  viscoplastic  part  A £vp  according  to: 

A£  =  AeE  +  Aevp  (10) 

The  stress  increment  is  related  to  the  elastic  strain  by  Hooke’s  law 

Aa  =  C  :  Ase  =  C  :  (Ae  -  A£vp)  (11) 

where  C  is  the  fourth-order  stiffness  tensor.  Next,  a  trial  elastic  stress:  dtriai  =  ot  +  C:Ae  is  cal¬ 
culated.  If / (dtriai,  svp,  T)  <  0,  the  increment  is  elastic  and  (jt+At  =  atriai-  If / (dtriai,  £Vp,  T)  >  0, 
then  there  is  viscoplastic  flow  and  the  stress  update  is  obtained  by  estimating  the  viscoplastic 
strain  increment  Asvp  (Simo  and  Hughes,  1998).  A  fully  implicit  integration  scheme  was  used, 
i.e.  the  viscoplastic  flow  was  determined  by  the  gradient  of  the  potential  function  at  t  +  At: 

Aevp  =  As'fAt  =  y<rA'>  A  t  (J £)  (12) 

Next,  Eq.  (12)  is  cast  into  the  following  form  (Wang  et  al.,  1997): 

a£vp  =  m(-)^  (13) 

where  the  multiplier  Ah  was  introduced  under  the  constraint  that  the  residual  Ar  given  by: 
Ah 

Ar  =  -£  -  (j)((7t+A\  ht+A\  T)  (14) 

tends  to  zero  upon  convergence  of  a  local  iterative  procedure.  Since  the  effective  stress  a  is 
homogeneous  of  degree  one  in  stresses,  A£vp  =  Ah.  Substituting  Eq.  (13)  into  Eq.  (11) 
yields: 
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Ag  =  C  : 


^Ae  -  Ah 


(15) 


To  compute  Ah,  a  Taylor  expansion  of  the  residual  is  considered  around  its  value  at  the 
previous  local  iteration,  i.e.: 


t+At  t+At 
'n+1  r n 


dr 


dAh 


t+At 


5h, 


n+1 


t+At 


:  Sg 


n+l 


dr 

dT 


t+\t 


ST, 


n+1 


(16) 


where  n  is  a  counter  for  the  local  iteration  (n  =  0  denotes  the  trial  elastic  state),  Shn+ \  = 
Ahn+\  —A hn,  Son+\  =  Aon+\  —  Aon,  and  5Tn+\  =  ATn+\  —  ATn.  Replacing  Eq.  (14)  in 
Eq.  (16)  we  obtain  (omitting  supra-indice): 


rn+ 1 


yAtShn+l 


dip 

da 


Sgh+ 1 


dcp 

dh 


Shn+ 1 


dcp 
d T 


st ; 


n+1 


(17) 


The  stress  update  for  the  trial  state  ( n  =  0)  is  found  by  assuming  that  the  total  strain 
increment  was  elastic,  while  for  n  >  0,  the  variation  of  Act  can  be  determined  by  consider¬ 
ing  the  effect  of  the  viscoplastic  strain  on  the  state  of  stress  as 


Sgu+i  =  —  C  :  Sh 


n+1 


n+1 


The  term  [||]  may  be  approximated  by: 


n+1 
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da 

dG 
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rd2  g1 
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:  Sg , 


n+1 
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ShS<T 


<5+} 

n  ) 


(18) 


(19) 


where  0  (0  ^  (9  <  1)  is  an  interpolation  parameter.  Combining  Eqs.  (18)  and  (19): 


^^n+l  H 


^  )  +  05hn 


8+ 

dhda 


Sh, 


n+1 


(20) 


where 


H  1 


=  c  1 


+  0Shn 


dG2 


(21) 


Further  substitution  of  Eq.  (20)  into  Eq.  (17)  yields 


Shn+ 1  p 


+ 


STn+\ 


with 


P 


+  0dhn 


a2q~ 

dhdG 


d(p  |  1 
dh  yAt 


(22) 


(23) 


In  Eq.  (22)  note  that,  since  the  total  strain  increment  is  fixed  during  the  local  iterative 
procedure,  Ssn+i  =  0.  The  derivatives  appearing  in  the  above  equations  are  calculated 
based  on  the  flow  rule  (see  Eq.  (2)): 
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(24a) 
(24b) 
(24c) 

4.  Results 

4.1.  Quasi-static  deformation  of  zirconium 
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In  this  section  and  the  following  we  present  applications  of  the  proposed  model  to  high- 
purity  zirconium.  Results  of  uniaxial  compression  tests  on  a  clock-rolled  zirconium  plate 
were  reported  by  Maudlin  et  al.  (1999b),  Kaschner  and  Gray  (2000),  and  Tome  et  al. 
(2001).  This  hep  aggregate  is  highly  anisotropic,  both  at  single  crystal  and  poly  crystal  lev¬ 
els.  It  was  processed  through  a  series  of  clock-rolling  and  annealing  cycles  to  produce  a 
plate  with  strong  basal  texture  «c)-axes  of  the  hep  crystals  predominantly  oriented  along 
the  plate  normal  direction).  The  process  of  multiple  rolling  passes  with  rotation  was  used 
in  order  to  obtain  a  nearly  isotropic  in-plane  texture.  Let  (x,y,z)  be  a  coordinate  reference 
frame  associated  with  the  plate,  where  the  (x,y)  plane  corresponds  to  the  original  rolling 
plane,  while  the  z-direction  is  the  through- thickness  (TT)  direction.  Right-circular  cylin¬ 
drical  tests  specimens  were  cut  from  the  plate  both  in  the  TT  (z— )  and  in-plane  (IP) 
(x-  or  y—)  directions.  Quasi-static  tests  in  both  IP  and  TT  directions  were  carried  out 
at  room  temperature  (298  K)  and  liquid  nitrogen  temperature  (77  K)  (Tome  et  al., 
2001).  In  this  paper  we  will  restrict  our  analysis  to  room  temperature  results.  It  was  found 
that  the  mechanical  response  is  highly  dependent  on  the  predominant  orientation  of  the 
(c)-axis  with  respect  to  the  loading  direction  (Tome  et  al.,  2001).  The  large  variations 
observed  in  flow  stress  and  strain-hardening  rate  are  due  to  the  relative  orientation  of 
the  available  twinning  and  slip  systems  and  the  critical  stress  levels  required  to  activate 
them.  In  the  case  of  TT  samples,  the  loading  is  applied  parallel  to  the  predominant  (c)-axis 
orientation.  Deformation  along  the  (c)-axis  requires  the  activation  of  {101 2}  (1011)  ten¬ 
sile  twinning,  {10ll}(TT23)  pyramidal  (c  +  a)-slip  or  the  much  harder  {112 2} (112 3) 
compressive  twinning.  Analysis  of  the  activity  of  these  deformation  systems  with  the 
VPSC  model  shows  that  at  room  temperature  tensile  twinning  and  pyramidal  (c  +  a)-slip 
dominate  in  tension  and  compression  along  the  (c)-axis,  respectively  (see  Fig.  6  in  Tome 
et  al.,  2001),  although  some  minor  activity  of  {1 1 2 2} ( 1 123)  compressive  twinning  and 
{1 121}(1 126)  tensile  twinning  has  been  also  observed  (Henrie  et  al.,  2004).  In  the  case 
of  IP  samples,  the  material  can  deform  within  the  plane  without  having  much  change  in 
dimension  along  the  harder-to-deform  (c)-axis.  Plastic  flow  occurs  at  lower  applied  stresses 
because  {10l0}(T2T0)  prismatic  slip  and  tensile  twinning  are  easily  activated  (see  also 
Fig.  6  of  Tome  et  al.  (2001),  which  shows  the  evolution  of  the  relative  activity  of  deforma¬ 
tion  modes  and  twin  fractions  with  deformation,  for  TT  and  IP  samples). 

In  what  follows,  all  the  reported  numerical  tests  using  the  VPSC  model  were  carried  out 
assuming:  (a)  the  experimental  initial  texture  (consisting  of  377  orientations)  of  the  clock- 
rolled  zirconium  plate,  (b)  the  deformation  mechanisms  characterized  operational  at  room 
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temperature  (i.e.  prismatic  (a)-slip,  pyramidal  (c  +  a)-slip  and  tensile  twinning),  and  (c)  the 
values  of  the  slip  and  twinning  parameters  (critical  stresses,  hardening  coefficients,  and 
rate  sensitivity  exponent),  as  reported  by  Tome  et  al.  (2001).  Comparison  between  the 
experimental  textures  and  textures  predicted  with  the  VPSC  model  were  provided  in  the 
same  paper. 

As  already  mentioned,  in  order  to  experimentally  characterize  the  anisotropy  induced  by 
texture  evolution,  the  material  should  be  subjected  to  a  specific  loading,  up  to  pre¬ 
determined  deformation  levels.  The  choice  of  the  pre-strain  loading  direction  should  be 
consistent  with  the  main  deformation  path  that  the  aggregate  will  undergo  in  a  specific 
application  (e.g.  uniaxial  compression  for  Taylor  impact  testing).  Then,  specimens  should 
be  cut  from  the  pre-strained  sample  and  further  tests  performed  to  quantify  the  anisotropy 
induced  by  deformation.  However,  this  would  require  a  large  number  of  mechanical  tests, 
which  could  be  very  difficult  to  carry  out.  Indeed,  there  is  little  data  for  large-strain  in-plane 
compression  of  Zr  in  sheet  form  because  measurements  are  limited  by  buckling.  Alterna¬ 
tively,  the  evolution  of  the  yield  surface  can  be  studied  by  performing  numerical  tests  using 
VPSC  model.  For  the  clock-rolled  plate  of  zirconium  described  previously,  the  VPSC  code 
was  first  used  to  characterize  the  evolution  of  the  yield  surface  during  IP  compression.  For 
this  purpose,  the  polycrystal  was  pre-strained  up  to  given  deformation  levels.  Fig.  1  shows 
the  biaxial  stress  plane  projections  of  the  VPSC  yield  loci  (symbols)  corresponding  to  pre¬ 
strain  levels  of  £*p  =  0.2%,  e^p  =  1%,  e^p  =  5%,  ^p  =  25%,  e^p  =  35%,  ^p  =  45%,  and 

=  60%.  Next,  for  each  individual  pre-strain  level  we  obtained  aJ  and  determined  the 
coefficients  L(e(p), ),a(&  )  (using  36  data  points  calculated  with  VPSC  and  a  least- 
square  fitting  procedure  in  conjunction  with  Eqs.  (5)— (7))  and  Yj  =  7(e(p)  as  well.  The 
homogeneity  parameter  was  considered  to  have  a  fixed  value  a  —  2  because  of  the  nearly 
elliptical  shape  of  the  yield  loci  for  high-purity  zirconium.  This  procedure  was  repeated 
for  all  seven  pre-strain  levels.  The  numerical  values  of  these  coefficients  for  each  CPB06 
yield  surface  are  listed  in  Table  1.  These  coefficients  quantify  the  very  strong  plastic  anisot¬ 
ropy  for  this  material  and  its  strong  evolution  with  strain.  Also,  the  linear  interpolation 
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Fig.  1.  Theoretical  yield  surface  evolution  for  a  zirconium  clock-rolled  plate  during  in-plane  compression  for 
various  levels  of  pre-strain.  Solid  lines:  CPB06  yield  surfaces  corresponding  to  fixed  pre-strain  levels,  determined 
using  the  VPSC  code  (symbols).  Dashed  lines:  CPB06  yield  surfaces  for  intermediate  pre-strain  levels  obtained  by 
linear  interpolation. 
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Table  1 


CPB06  yield  function  coefficients  corresponding  to  various  pre-strain  levels  for  uniaxial  compression  of  clock- 
rolled  zirconium  along  an  in-plane  direction 


k 

L\2 

L\3 

T22 

T23 

£33 

L44 

L55 

^66 

0.2% 

0.071 

2.883 

1.875 

1.252 

1.965 

1.204 

0.989 

1.075 

3.956 

1% 

0.224 

2.15 

1.407 

1.179 

1.438 

0.892 

0.604 

0.617 

1.691 

5% 

0.500 

2.923 

1.823 

0.865 

1.745 

0.418 

0.984 

1.145 

2.980 

25% 

0.263 

2.704 

1.727 

1.821 

2.056 

0.869 

0.949 

0.791 

1.403 

35% 

0.174 

2.505 

1.683 

2.320 

2.054 

1.012 

0.889 

0.694 

1.098 

45% 

-0.269 

0.853 

1.444 

1.143 

0.624 

1.007 

0.638 

0.461 

0.665 

60% 

-0.220 

0.766 

0.858 

0.637 

0.900 

0.644 

0.299 

0.189 

0.255 

For  all  levels  of  effective  plastic  strain:  a  =  2  and  Ln  =  1.0 


procedure  was  used  to  obtain  the  macroscopic  yield  surfaces  corresponding  to  any  pre¬ 
strain  level  (see  Eq.  (8)).  In  Fig.  1,  the  dashed  lines  represent  the  CPB06  yield  surfaces 
for  intermediate  pre-strain  levels  obtained  using  linear  interpolation.  Note  that  an  accurate 
description  of  the  evolving  yield  function  for  pre-strain  levels  of  around  30%  and,  further, 
up  to  60%  are  needed  for  the  quasi-static  IP  compression  simulation  below,  and  for  the 
Taylor  cylinder  simulation  of  next  subsection,  respectively.  It  is  also  worth  mentioning  that 
the  successive  yield  function  projections  for  pre-strains  higher  than  5%  are  more  rounded 
and  centered,  reflecting  a  decrease  of  the  plastic  anisotropy  (note  that  all  anisotropy  coef¬ 
ficients  have  initially  rather  large  values  which  decrease  at  large  strains).  This  effect  is 
related  to  the  high  tensile  twinning  activity  and  the  concomitant  strong  reorientation  of 
the  twinned  portions  of  the  grains  that  take  place  during  in-plane  compression. 

The  proposed  model  was  used  to  simulate  the  quasi-static  uniaxial  compression  of  a 
right  cylinder  of  circular  cross  section  with  its  axis  parallel  to  the  in-plane  x-direction 
up  to  28%  longitudinal  strain,  for  comparison  with  analog  experimental  results  reported 
by  Tome  et  al.  (2001)  (Fig.  2).  Note  that  because  the  z-direction  is  very  hard  to  deform, 


Fig.  2.  Cross-sections  for  quasi-static  compression  of  a  zirconium  cylinder  with  axis  parallel  to  an  in-plane 
direction,  (a)  Final  ABAQUS  FE  mesh  obtained  with  the  proposed  model,  (b)  Photograph  of  the  final 
experimental  cross-section  and  EPIC/VPSC  simulation  results  (dashed  lines),  after  Tome  et  al.  (2001).  The  arrows 
indicate  the  predominant  direction  of  the  crystallographic  (c)-axes. 
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the  specimen  deforms  mostly  in  the  y-direction  and  the  cross-section  of  the  deformed  spec¬ 
imen  becomes  elliptical.  Our  simulations  predict  in-plane  expansion  of  22%  and  out-of- 
plane  expansion  of  6%  (or  an  ovalization  long/short  radii  =  1.19),  which  agrees  very  well 
with  the  experimental  results  reported  in  Tome  et  al.  (2001).  It  is  worth  noting  that  under 
the  assumption  of  isotropic  hardening  and  for  fixed  values  of  the  anisotropy  coefficients 
(i.e.  texture  evolution  was  neglected  and  the  anisotropy  coefficients  for  any  pre-strain  level 
were  set  equal  to  the  values  of  L  obtained  based  on  VPSC  yield  stresses  calculated  by 
probing  the  0.2%  pre-strained  polycrystal  in  different  directions),  the  predicted  expansions 
would  have  been  of  27%  in-plane  and  1%  out-of-plane.  This  clearly  implies  that  a  good 
agreement  with  the  experimental  results  can  be  obtained  only  if  the  changes  in  anisotropy 
induced  by  texture  evolution  (in  this  case,  a  decrease  in  plastic  anisotropy)  is  taken  into 
account.  Fig.  2a  shows  the  final  ABAQUS  FE  mesh.  The  photograph  of  the  final  exper¬ 
imental  cross-section  shape  and  the  superimposed  simulation  result  (obtained  using  the  FE 
program  EPIC  directly  linked  to  the  VPSC  code)  are  shown  in  Fig.  2b  (simulation  and 
data  after  Tome  et  al.,  2001).  The  comparison  demonstrates  that  very  good  agreement 
with  the  experimental  results  can  be  obtained  by  using  an  appropriate  initial  yield  surface 
representation  (CPB06)  and  a  hardening  law  that  accounts  for  texture  evolution. 

Next,  we  have  simulated  the  quasi-static  through-thickness  compression  of  the  zirco¬ 
nium  specimen.  Fig.  3  shows  the  yield  stresses  obtained  with  the  VPSC  code  (symbols) 
for  five  individual  levels  of  compressive  pre-strain  (0.2%,  1%,  5%,  25%  and  35%)  along 
the  through-thickness  direction,  the  corresponding  five  individual  CPB06  yield  surfaces 
for  which  the  coefficients  were  determined  (solid  lines),  and  yield  surfaces  at  intermediate 
levels  of  strain  obtained  by  linear  interpolation  (dashed  lines).  Note  that  a  description  of 
the  evolving  yield  function  up  to  pre-strain  levels  of  around  35%  was  enough  for  the  28% 
strain  quasi-static  TT  compression  simulation  below. 


Fig.  3.  Theoretical  yield  surface  evolution  for  a  zirconium  clock-rolled  plate  during  through- thickness 
compression  for  various  levels  of  pre-strain.  Solid  lines:  CPB06  yield  surfaces  corresponding  to  fixed  pre-strain 
levels,  determined  using  the  VPSC  code  (symbols).  Dashed  lines:  CPB06  yield  surfaces  for  intermediate  pre-strain 
levels  obtained  by  linear  interpolation. 
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Table  2 


CPB06  yield  function  coefficients  corresponding  to  various  pre-strain  levels  for  uniaxial  compression  of  clock- 
rolled  zirconium  along  the  through-thickness  direction 


k 

L\2 

L\2> 

T22 

T23 

£33 

L44 

L55 

^66 

0.2% 

0.07 

2.883 

1.875 

1.252 

1.965 

1.204 

0.989 

1.075 

3.956 

1% 

0.23 

2.964 

1.556 

0.803 

1.476 

0.433 

1.152 

1.173 

3.313 

5% 

0.40 

2.652 

1.619 

0.813 

1.541 

0.583 

0.988 

1.127 

2.858 

25% 

0.19 

3.348 

1.575 

1.039 

1.525 

0.603 

1.883 

2.081 

3.391 

35% 

0.08 

3.272 

1.524 

1.120 

1.471 

0.641 

2.122 

2.312 

3.847 

For  all  levels  of  effective  plastic  strain:  a  =  2  and  Ln  =  1.0. 


The  numerical  values  of  the  anisotropy  coefficients  for  CPB06  corresponding  to  each 
individual  level  of  pre-strain  are  listed  in  Table  2.  Note  that  the  CPB06  yield  surface  evo¬ 
lution  associated  to  deformation  induced  anisotropy  in  uniaxial  compression  is  very  differ¬ 
ent  than  the  evolution  associated  to  uniaxial  compression  in  the  in-plane  direction  (cf. 
Figs.  1  and  3).  Through- thickness  compression  is  dominated  by  slip.  Slip-generated  crystal 
rotation  would  normally  result  in  texture  evolution  at  significant  levels  of  strain.  However, 
in  this  case  the  texture  is  not  expected  to  evolve  much  from  the  initial  stable  strong  basal 
texture.  On  the  other  hand,  during  in-plane  compression  tensile  twinning  is  easily  activated 
thus  producing  a  strong  evolution  in  texture. 

The  yield  surface  evolution  depicted  in  Fig.  3  was  in  turn  used  in  an  ABAQUS  FE  sim¬ 
ulation  of  compression  of  a  right  cylinder  of  circular  cross  section  and  axis  parallel  to  the 
z-axis  up  to  28%  longitudinal  strain,  and  compared  with  the  experimental  results  (Tome 
et  al.,  2001)  (Fig.  4).  Since  compression  takes  place  along  the  axis  of  symmetry  of  the  sam¬ 
ple,  the  final  cross-section  remains  circular.  The  proposed  model  correctly  predicts  this 
experimental  trend  (cf.  Fig.  4a  and  b). 

4.2.  High  strain-rate  deformation  of  zirconium 

Controlled  impact  tests  are  used  to  characterize  the  high  strain-rate  mechanical 
response  of  metals.  Such  tests  provide  results  in  the  104— 105  s-1  strain-rate  range,  directly 


Fig.  4.  Cross-sections  for  quasi-static  compression  of  a  zirconium  cylinder  with  axis  parallel  to  the  through¬ 
thickness  direction,  (a)  Final  ABAQUS  FE  mesh  obtained  with  the  proposed  model,  (b)  Photograph  of  the  final 
experimental  cross-section  and  EPIC/VPSC  simulation  results  (dashed  lines),  after  Tome  et  al.  (2001). 
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above  that  data  accessible  via  Hopkinson  Bar  techniques  and  below  shock  experimenta¬ 
tion.  One  standard  approach  used  is  the  impact  of  a  right  circular  cylinder  on  either  a  mas¬ 
sive  anvil  (rigid  plate)  or  a  receptor  rod  of  identical  material.  This  experiment  is  commonly 
called  a  Taylor  impact  test  (Taylor,  1948).  Due  to  the  impact,  the  rod  plastically  deforms 
and  shortens  causing  the  material  at  the  impact  surface  to  flow  radially  outward  relative  to 
the  rod  axis. 

Maudlin  et  al.  (1999b)  reported  results  of  Taylor  cylinder  impact  tests  using  in-plane 
specimens  cut  from  the  same  zirconium  plate  described  and  characterized  in  the  previous 
section  (i.e.  cylindrical  axes  were  coincident  with  in-plane  directions).  The  specimens  were 
50.8  mm  in  length  and  have  a  length  over  diameter  ratio  of  6.67.  The  anvil  was  AF1410 
steel,  heat-treated  to  a  surface  hardness  of  Rc  58  and  lapped  to  a  mirror  finish.  The  cylin¬ 
ders  were  fired  using  a  gas-driven  gun  at  various  velocities:  50,  101,  170,  and  243  m/s.  Only 
post- test  data  for  the  243  m/s  impact  velocity  were  reported.  As  discussed  in  Maudlin  et  al. 
(1999b),  the  data  set  consists  of  high-resolution  scans  of  the  photographed  footprint,  i.e. 
the  cross-sectional  area  at  the  impact  interface  (on  the  y-z-plane)  and  of  the  photographed 
side  profiles  for  the  major  and  minor  dimensions.  From  the  analysis  of  the  digitized  side 
profiles,  axial  strain  profiles  (log  strain)  as  a  function  of  axial  position,  measured  relative 
to  the  impact  interface,  were  obtained. 

We  have  applied  the  proposed  elastic-viscoplastic  model  to  the  description  of  the  high 
strain-rate  deformation  of  the  clock-rolled  zirconium  plate  at  room  temperature  and  com¬ 
pared  the  simulations  with  Taylor  impact  test  results  reported  by  Maudlin  et  al.  (1999b). 
In  addition  to  the  yield  surface  evolution  due  to  plastic  deformation,  the  viscosity  coeffi¬ 
cient  y  and  the  strain-rate  sensitivity  parameter  m  (see  Eqs.  (2)  and  (3))  needed  to  be  deter¬ 
mined.  Results  of  dynamic  in-plane  tests  performed  on  the  clock-rolled  Zr  material  at 
strain-rates  from  1000/s  to  3500/s  conducted  using  a  Split  Hopkinson  Pressure  Bar 
(SHPB)  were  reported  by  Kaschner  and  Gray  (2000).  Based  on  these  tests,  we  have  deter¬ 
mined  y  =  2500  s-1  and  m  =  7.0  (see  Fig.  5).  The  constitutive  behavior  for  the  in-plane 
compression  corresponding  to  strains  beyond  the  range  of  the  experimental  data  was  esti¬ 
mated  based  on  polycry stalline  calculations. 


Fig.  5.  Stress-strain  curves  of  in-plane  compression  of  clock-rolled  zirconium  samples  for  various  strain-rates. 
Symbols:  experimental  data,  after  Kaschner  and  Gray  (2000).  Lines:  simulation  results  using  the  proposed  elastic- 
viscoplastic  model  for  y  =  2500  s-1  and  m  =  7.0. 


B.  Plunkett  et  al.  /  International  Journal  of  Plasticity  23  (2007)  1001-1021  1013 

In  order  to  simulate  Taylor  impact  tests  of  zirconium,  the  proposed  elastic- viscoplastic 
model  was  implemented  into  an  ABAQUS  explicit  user  material  subroutine.  The  rate- 
independent  yield  surface  characterization  was  the  one  described  in  the  previous  subsec¬ 
tion,  based  on  VPSC  calculations  of  the  induced  anisotropy  during  quasi-static  in-plane 
compression  coupled  with  CPB06  (see  Fig.  3).  It  should  be  noted  that  by  making  this 
choice,  we  have  disregarded  the  possibility  of  changes  in  the  combination  of  active  defor¬ 
mation  mechanisms  at  single  crystal  level  as  strain-rate  increases.  For  comparison  pur¬ 
poses,  simulations  under  the  assumption  of  isotropic  hardening  and  anisotropy 
coefficients  set  equal  to  the  anisotropy  coefficients  corresponding  to  0.2%  pre-strain  (tex¬ 
ture  evolution  neglected)  were  also  carried  out. 

The  zirconium  cylinder  was  modeled  using  2117  ABAQUS  C3D8R  linear  brick  ele¬ 
ments  with  free  boundary  conditions.  The  anvil  target  was  modeled  as  an  analytical  rigid 
surface.  Simulations  were  done  for  the  impact  velocity  of  243  m/s.  The  step  size  used  to 
simulate  the  event  was  At  =  2  x  10-8  s.  Calculations  were  performed  for  90  (is,  at  which 
time  the  specimen  had  rebound  off  of  the  target  and  all  plastic  deformation  had  ceased. 
Fig.  6  shows  a  comparison  between  the  experimental  and  simulated  logarithmic  strain  pro¬ 
file  along  the  major  and  minor  axes  of  the  post-test  specimen.  The  strains  along  the  major 
and  minor  elliptical  axes,  which  are  both  zero  at  large  z,  accumulate  at  different  rates  as  z 
decreases.  A  good  agreement  with  the  experimental  data  is  observed.  Fig.  7  depicts  a  visual 
comparison  between  the  simulated  and  experimental  major  and  minor  profiles  and  foot- 


Fig.  6.  Comparison  of  theoretical  (lines)  and  experimental  (symbols)  logarithmic  strain  profiles  of  the  post-test 
zirconium  Taylor  impact  specimen.  Data  after  Maudlin  et  al.  (1999b). 


1014 


B.  Plunkett  et  al.  /  International  Journal  of  Plasticity  23  (2007)  1001-1021 


Fig.  7.  Comparison  of  the  simulated  (with  contours  of  axial  plastic  strain)  and  experimental  cross-sections  of  the 
post-test  zirconium  Taylor  impact  experiment  for  the  major  profile,  the  minor  profile,  and  the  footprint.  Data 
after  Maudlin  et  al.  (1999b). 


print  of  the  post-test  specimen.  The  slight  differences  between  the  experimental  data  and 
the  simulation  results  are  most  likely  related  to  temperature  effects  which  were  not  taken 
into  consideration  in  the  simulations  because  experimental  data  was  not  available.  Note 
that  the  geometry  of  the  footprint  of  the  specimen  clearly  indicates  a  decrease  in  anisot¬ 
ropy  (the  ratio  of  major  to  minor  diameter  is  «1.1).  This  is  due  to  the  large  twinned  vol¬ 
ume  fraction  resulting  from  strains  of  around  60-70%  at  the  impact  section  (see  strain 
levels  in  Fig.  7)  which  tends  to  decrease  the  intensity  of  the  texture  (as  already  seen  in 
Fig.  1  and  discussed  in  Section  4.1). 

As  discussed  earlier  in  the  quasi-static  case,  the  consideration  of  a  hardening  law  that 
accounts  for  texture  evolution  is  essential  to  obtain  the  good  agreement  with  the  experi¬ 
mental  results.  Fig.  8  shows  a  comparison  between  experimental  minor  and  major  profiles 
of  the  post-test  specimen  and  those  calculated  using  isotropic  hardening  and  fixed  values 
for  the  anisotropy  coefficients  corresponding  to  0.2%  pre-strain.  It  is  evident  that  the 
effects  of  substantial  twinning  activity  can  be  captured  only  by  a  hardening  law  that 
reflects  the  evolution  of  texture. 

4.3.  High  strain-rate  deformation  of  tantalum 

The  anisotropic  elastic-viscoplastic  model  was  originally  developed  to  capture  the 
anisotropic,  asymmetric  hardening  behavior  of  hep  metals.  The  static  yield  surface  in  this 
model  is  the  CPB06  yield  surface.  However,  the  CPB06  criterion  can  be  also  used  to  rep- 
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Fig.  8.  Comparison  of  theoretical  (lines)  and  experimental  (symbols)  logarithmic  strain  profiles  of  the  post-test 
zirconium  Taylor  impact  specimen.  Simulation  were  performed  with  an  isotropic  hardening  law.  Data  after 
Maudlin  et  al.  (1999b). 

resent  the  yielding  behavior  of  a  material  with  cubic  crystallographic  structure  deforming 
solely  by  slip  by  simply  setting  the  strength-differential  coefficient  k  to  zero.  In  this  section, 
we  apply  the  proposed  formulation  to  the  description  of  the  anisotropy  of  a  bcc  tantalum 
sheet  (obtained  by  unidirectional  rolling  and  displaying  a  near-orthotropic  texture)  sub¬ 
jected  to  high  strain-rate  and  large  plastic  deformations  (data  after  Maudlin  et  al., 
1999a, b). 

Unlike  for  the  previously  studied  clock-rolled  zirconium,  in  this  case  there  was  no  suf¬ 
ficient  experimental  data  available  to  quantify  the  deformation  induced  anisotropy  using 
VPSC  simulations  (i.e.  no  data  based  on  which  we  could  determine  the  parameters 
involved  in  the  VPSC  model  were  available).  Therefore,  to  estimate  the  coefficients 
involved  in  the  macroscopic  CPB06  yield  criterion,  we  have  used  the  values  of  yield  stres¬ 
ses  in  different  directions  reported  by  Maudlin  et  al.  (1999a,b).  The  numerical  values  of  the 
CPB06  coefficients  are  listed  in  Table  3.  Moreover,  due  to  lack  of  information  regarding 
the  yield  behavior  in  pure  shear,  it  was  assumed  a  value  equal  to  0.577  times  the  yield 
stress  in  uniaxial  tension  along  the  rolling  direction.  Fig.  9  shows  the  projection  in  the 
biaxial  plane  of  the  CPB06  orthotropic  yield  surface  for  tantalum  (solid  lines)  together 
with  the  yield  stresses  (symbols)  calculated  by  Maudlin  et  al.  (1999a,b)  using  a  Taylor- 
Bishop-Hill  (TBH)  polycrystal  model. 

Uniaxial  compression  test  results  at  various  strain-rates  and  temperatures  were 
reported  for  tantalum  specimens  cut  from  the  same  rolled  plate.  Results  of  Split-Hopkin- 
son  bar  tests  performed  at  different  temperatures  in  the  range  200-1000  °C  and  for  strain- 
rates  in  the  range  2800-3900  s-1  have  shown  that  the  compressive  response  depends  on  the 
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Table  3 


CPB06  yield  function  coefficients  corresponding  to  rolled  tantalum;  a  =  2  and  Ln  =  1.0 


k 

Ti2 

Lu 

Lll 

L23 

L33 

0 

-0.1911 

-0.0687 

1.0411 

0.0067 

1.1366 

Fig.  9.  Line:  CPB06  yield  surface  corresponding  to  a  rolled  sheet  of  tantalum.  Symbols:  yield  stress  points 
predicted  with  the  Taylor-Bishop-Hill  approach.  Stresses  are  normalized  about  the  yield  stress  in  the  x-direction 
(rolling  direction).  Data  after  Maudlin  et  al.  (1999a,b). 


strain-rate  and  temperature.  The  yield  stress  was  found  to  increase  with  increasing  strain- 
rate  and  decreasing  temperature.  To  capture  the  effect  of  temperature  on  the  compressive 
yield  stress  of  tantalum  we  used  a  Johnson-Cook  type  hardening  law  (Johnson  and  Cook, 
1983): 


r(fivp,  t) 


To (^vp) 


T~Tvm 
Tmelt  -  Tn 


(25) 


In  Eq.  (25),  T0(eVp)  is  the  rate-independent  hardening  law  at  room  temperature,  Trm  is 
the  room  temperature  (typically  25  °C),  Tme\t  is  the  melting  temperature,  and  tm  is  a  mate¬ 
rial  constant.  For  tantalum,  Tmeit  =  2977  °C,  while  a  value  of  0.42  for  the  exponent  tm 
gives  an  overall  good  approximation  of  the  experimental  data.  Assuming  that  all  mechan¬ 
ical  work  is  converted  into  heat,  from  the  energy  equation  (Malvern,  1969)  the  tempera¬ 
ture  rise  due  to  mechanical  deformation  is: 


A  Tm 


g  :  A8vp 

pCp 


(26) 


For  tantalum,  p  =  16,640  kg/m3,  <^=140J/kgK.  Furthermore,  using  the  available 
data,  a  viscosity  y  =  1500  s-1  and  a  strain-rate  sensitivity  parameter  m  =  5.0  (see  Eq. 
(2))  were  determined. 
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Next,  we  evaluated  the  predictive  capability  of  the  present  model  to  simulate  the 
severe  geometry  changes  of  the  tantalum  Taylor  specimens  reported  by  Maudlin  et  al. 
( 1999a, b).  The  cylinders  were  38.1  mm  long  and  7.6  mm  in  diameter.  The  longitudinal 
axis  of  the  specimen  was  either  parallel  to  the  rolling  direction  (y)  or  transverse  direction 
(x).  Only  results  on  specimens  having  the  longitudinal  axis  parallel  to  the  rolling  direc¬ 
tion  were  reported.  The  cylinders  were  launched  using  a  powder  gun.  The  velocity  of 
the  cylinders  measured  using  both  pressure  transducers  and  parallel  laser  beams  were 
175  m/s  ±3  m/s.  The  anvil  target  was  made  of  steel.  Metallographic  analysis  of  the 
pre-  and  post-test  specimens  showed  negligible  texture  evolution.  Therefore,  the  experi¬ 
ment  was  modeled  assuming  isotropic  hardening.  Because  the  material  was  orthotropic, 
only  a  quarter  section  was  modeled  using  8170  ABAQUS  C3D8R  linear  brick  elements 
with  free  boundary  conditions  and  an  initial  velocity  of  175  m/s.  The  step  size  used  to 
simulate  the  event  was  At  =  2x  10-8  s  and  was  conducted  until  100  jlls.  At  100  jlls  the 
specimen  had  rebound  off  of  the  target,  and  all  plastic  deformation  had  ceased.  The  anvil 
target  was  modeled  as  an  analytical  rigid  surface.  A  higher  number  of  elements  were  used 
for  the  tantalum  simulations  as  compared  to  the  zirconium  simulations  because  the  tan¬ 
talum  specimen  undergoes  larger  deformations  (for  tantalum  at  z  =  0,  the  strain  along 
the  major  and  minor  dimensions  were  of  87%  and  70%,  respectively  while  for  zirconium 
the  values  of  the  strains  at  z  =  0  were  of  44%  and  31%,  respectively).  Axial  strain  profiles 
for  both  the  major  and  minor  sides  of  the  deformed  specimens  (cylinder  axis  parallel  to 
the  rolling  direction)  in  comparison  with  experiments  are  shown  in  Fig.  10.  Very  good 
agreement  is  observed.  The  simulated  and  experimental  footprints  of  the  post  test  spec¬ 
imen  are  shown  in  Fig.  11.  The  data  shown  in  Figs.  10  and  11  are  the  results  from  three 
separate  experiments  conducted  for  different  specimens  under  the  same  conditions.  The 
simulated  and  experimental  post- test  specimens  are  visually  compared  in  Fig.  12.  In 


Fig.  10.  (a)  Major  and  (b)  minor  profiles  of  the  tantalum  Taylor  impact  specimen.  Solid  lines  represent 
simulation  results,  symbols  represent  experimental  data.  Data  after  Maudlin  et  al.  ( 1999a, b). 
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Fig.  11.  Footprint  of  the  tantalum  Taylor  impact  specimen.  Solid  lines  represent  simulation  results,  symbols 
represent  experimental  data.  Data  after  Maudlin  et  al.  (1999a,b). 
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Fig.  12.  Comparison  between  simulated  (with  contours  of  axial  plastic  strain)  and  experimental  tantalum  Taylor 
impact  specimens.  Photographs  of  the  post-test  specimen  after  Maudlin  et  al.  ( 1999a, b). 
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Maudlin  et  al.  (1999b),  the  same  tests  were  simulated  by  introducing  anisotropy  through 
direct  link  between  the  Taylor-Bishop-Hill  polycrystal  model  and  the  FE  code  EPIC, 
assuming  the  presence  of  a  polycrystal  at  each  integration  point.  Their  results  were  also 
in  very  good  agreement  with  the  experimental  data.  However,  Figs.  10-12  clearly  dem¬ 
onstrate  that  very  accurate  results  can  also  be  obtained  with  the  present  formulation, 
at  a  fraction  of  the  computation  time. 

5.  Conclusions 

A  macroscopic  anisotropic  elastic-viscoplastic  model  that  captures  the  influence  of 
evolving  texture  on  the  mechanical  response  of  textured  metals  for  both  quasi-static 
and  dynamic  loading  conditions  was  proposed.  Yielding  was  described  using  the  aniso¬ 
tropic  yield  function  proposed  by  Cazacu  et  al.  (2006).  This  yield  function  can  capture 
simultaneously  both  anisotropy  and  tension/compression  asymmetry.  Moreover,  if  the 
parameter  associated  with  strength  differential  effects  is  set  to  zero,  the  yield  function  also 
describes  with  accuracy  yielding  of  metals  with  cubic  symmetry. 

The  anisotropy  coefficients  as  well  as  the  size  of  the  elastic  domain  were  considered  to 
be  functions  of  the  accumulated  plastic  strain.  To  capture  the  anisotropy  induced  by  tex¬ 
ture  evolution,  a  multi-scale  methodology  consisting  in  using  experimental  textures  and 
stress-strain  curves,  VPSC  polycrystalline  calculations,  and  macroscopic  scale  interpola¬ 
tion  techniques  was  used.  The  procedure  started  by  fitting  the  parameters  of  the  poly  crys¬ 
tal  model  to  reproduce  the  mechanical  response  of  the  pre-textured  material.  Next,  using 
this  adjusted  polycrystal  model,  the  evolution  of  texture  was  calculated  for  fixed  pre-strain 
level  along  a  given  deformation  path.  Then,  in  order  to  quantify  the  induced  anisotropy, 
the  pre-strained  polycrystal  was  numerically  probed  along  various  directions.  The  calcu¬ 
lated  yield  stresses  along  these  directions  were  in  turn  used  to  determine  the  CPB06  coef¬ 
ficients.  This  procedure  was  repeated  for  a  finite  set  of  pre-strain  levels  and  linear 
interpolation  was  in  turn  used  to  obtain  the  macroscopic  yield  surfaces  corresponding 
to  an  arbitrary  pre-strain  level. 

The  above  procedure  was  used  to  model  the  quasi-static  mechanical  response  of  high- 
purity  zirconium.  Since  the  proposed  model  accounts  for  the  deformation  mechanisms 
associated  with  a  given  deformation  path,  it  captured  the  strong  anisotropy  and  texture 
evolution  produced  by  twinning  which  is  observed  when  the  material  is  subjected  to  in¬ 
plane  compression.  On  the  other  hand,  the  resulting  macroscopic  yield  surface  evolution 
along  through-thickness  compression  was  very  different  than  the  one  associated  with 
in-plane  compression  (cf.  Figs.  1  and  3).  Through-thickness  compression  is  dominated 
by  slip  so  the  texture  does  not  significantly  change  during  loading,  whereas  during 
in-plane  compression  tensile  twinning  is  easily  activated,  producing  a  strong  texture 
evolution. 

Based  on  the  above  model  for  quasi-static  behavior,  an  overstress  type  of  approach  was 
in  turn  used  to  account  for  rate  effects.  The  resulting  elastic-viscoplastic  model  was  applied 
to  the  description  of  the  high  strain-rate  response  of  low-symmetry  (clock-rolled  hep  zir¬ 
conium)  and  high-symmetry  (rolled  bee  tantalum)  textured  metals.  The  very  good  agree¬ 
ment  between  the  simulated  and  experimental  post-test  geometries  of  the  Taylor-impact 
specimens  in  terms  of  major  and  minor  side  profiles  and  impact  footprints  proved  the  abil¬ 
ity  of  the  model  to  describe  the  striking  difference  in  the  evolution  of  anisotropy  between 
zirconium  and  tantalum. 
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